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Abstract. We analyzed intermittency in the solar wind, as observed on the ecliptic plane, 
looking at magnetic field and velocity fiuctuations between 0.3 and 1 AU, for both fast 
and slow wind and for compressive and directional fluctuations. Our analysis focused on 
the property that probability distribution functions of a fluctuating field affected by in- 
termittency become more and more peaked at smaller and smaller scales. Since the peaked- 
ness of a distribution is measured by its flatness factor we studied the behavior of this 
parameter for different scales to estimate the degree of intermittency of our time series. 
We confirmed that both magnetic field and velocity fluctuations are rather intermittent 
and that compressive magnetic fluctuations are generally more intermittent than the cor- 
responding velocity fluctuations. In addition, we observed that compressive fiuctuations 
are always more intermittent than directional fluctuations and that while slow wind in- 
termittency does not depend on the radial distance from the sun, fast wind intermittency 
of both magnetic field and velocity fiuctuations clearly increases with the heliocentric 
distance. 

We propose that the observed radial dependence can be understood if we imagine in- 
terplanetary fiuctuations made of two main components: one represented by coherent, 
non propagating structures convected by the wind and, the other one made of propa- 
gating, stochastic fiuctuations, namely Alfven waves. While the first component tends 
to increase the intermittency level because of its coherent nature, the second one tends 
to decrease it because of its stochastic nature. As the wind expands, the Alfvenic con- 
tribution is depleted because of turbulent evolution and, consequently, the underlying 
coherent structures convected by the wind, strengthen further on by stream-stream dy- 
namical interaction, assume a more important role increasing intermittency, as observed. 
Obviously, slow wind doesn't show a similar behavior because Alfvenic fiuctuations have 
a less dominant role than within fast wind and the Alfvenicity of the wind has already 
been frozen by the time we observe it at 0.3 AU. Finally, our analysis suggests that the 
most intermittent magnetic fiuctuations are distributed along the local interplanetary 
magnetic field spiral direction while, those relative to wind velocity seem to be located 
along the radial direction. 



1. 1. Introduction 

The basic view that we have of the solar wind is that 
of a magnetofluid pervaded by fluctuations over a wide 
range of scales which are strongly modified by the effects of 
the dynamics during the expansion into the interplanetary 
medium. These effects are more relevant within the inner 
heliosphere and on the Ecliptic where the stream-stream 
dynamics more strongly reprocesses the original plasma and 
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the large velocity shears add new fluctuations to the origi- 
nal spectrum (Coleman, 1968, Roberts et al., 1991). This 
scenario has reconciled the "wave" point of view proposed 
by Belcher and Davis (1971), i.e. solar origin of the fluctu- 
ations, and the "turbulence" point of view, i.e. local gen- 
eration due to velocity shears, proposed by Coleman(1968). 
The first consequence of this scenario is that large fiuctua- 
tions of solar origin containing energy interact non-linearly 
with other fiuctuations of local origin giving rise to an en- 
ergy exchange between different scales, which can be inter- 
preted as the usual energy cascade towards smaller scales 
in fully developed turbulence. As a matter of fact, space- 
craft observations have shown that the spectral slope of the 
power spectrum of these fiuctuations changes with the ra- 
dial distance from the sun (Bavassano et., 1982, Denskat and 
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Neubaucr, 1983). This behaviour was recognized (Tu et al., 
1984) as a clear experimental evidence that cascade pro- 
cesses due to non-linear interaction between opposite prop- 
agating Alfven waves were active in the solar wind with. 
One of the consequences of this radial evolution was the 
observed radial decrease of the correlation of velocity and 
magnetic field fluctuations (generally known as cross helic- 
ity, or Alfvenicity) (Roberts et al. 1987). These observa- 
tions finally answered to the question of whether the ob- 
served fluctuations wore remnants of coronal processes or 
were dynamically created during the expansion. However, 
successive theoretical models (Sec review by Tu and Marsch, 
1995) which tried to obtain the radial spectral evolution of 
the solar wind fluctuations had to deal with peculiarities of 
the observations that they could not reproduce within the 
framework of solely non-linear interacting waves. The lack 
of a strict self-similarity of the fluctuations and the conse- 
quent non applicability of strict scale invariance (Marsch 
and Liu, 1993), the strong anisotropy shown by velocity 
and field fluctuations (Bavassano et al., 1982, Tu et al., 
1989, Roberts, 1992), the difl^erent radial evolution of the 
minirrmm variance direction for magnetic field and veloc- 
ity (Klein et al., 1993), the lack of equipartition between 
magnetic and velocity fluctuations (Matthaeus and Gold- 
stein, 1982, Bruno et al., 1985) all contributed to suggest 
the idea that fluctuations could possibly be due to a mix- 
ture of propagating waves and static structures convected 
by the wind. Some kind of filamentary structure, simi- 
lar to fiux tubes, was firstly proposed by McCracken and 
Ness (1966) and the observed spectral radial evolution of 
the large scale fluctuations has been attributed to the in- 
teraction of outward propagating Alfvcn waves with these 
structures (Tu and Marsch, 1993, Bruno and Bavassano, 
1991; Bavassano and Bruno, 1992). Incompressible mag- 
netic structures were found by Tu and Maxsch (1991)and 
magnetic fluctuations with a large correlation length par- 
allel to the ambient magnetic field, suggested the idea of a 
quasi-two-dimensional, incompressible turbulence for which 
k ■ B = (Matthaeus et al., 1990). Thus, solar wind fiuc- 
tuations are not isotropic and scale-invariant, two of the 
fundamental hypotheses at the basis of K41 Kolmogorov's 
theory (1941). This theory is based on an important statis- 
tical relation, which characterizes turbulent flows, between 
velocity increments 5vr =< \V{x + r) — V{x)\ >, measured 
along the flow direction x, and the energy transfer rate e 
at the scale separation r = |r|, that is 6vr ~ (er)^''^ or, 
more in general, Sv^ ~ {ery^^. If e is constant, the previous 
relation simply reads 6v? ^ r^^^ and fluctuations are said 
to be self-similar, and our signal is a simple fractal. How- 
ever, as remarked by Landau (Kolmogorov, 1962, Obukhov, 
1962), if e statistically depends on scale due to the mecha- 
nism that transfers energy from larger to smaller eddies, e 
will be replaced by er and a new scaling has to be evaluated 
5vP ~ rf/^ < >. Expressing e?"^ via a scaling rela- 
tion with r, we obtain < e?!^'* >~ r'^p/^ and, consequently, 
Sv^! ~ r"'' where Sp = p/3 + rp/3 is generally a nonlinear 
function of p. This means that the global scale invariance 
required in the K41 theory would release towards a local 
scale invariance where different fractal sets characterized by 
different scaling exponents can be found. 

One of the consequences of this la.ck of a universal scale 
invariance, directly observed in experimental tests, is that 
the shape of the probability density functions (PDFs) of the 
velocity increments at a given scale is not the same for each 
scale but roughly evolves from a Gaussian shape, near the 



integral scale, to a distribution whose tails are much flatter 
than those of a Gaussian, resembling a stretched exponen- 
tial near the dissipation scale. This means that the largest 
events, contained in the tails of the distribution, do not fol- 
low the Gaussian statistics but show a much larger prob- 
ability. This phenomenon is also called intermittency and, 
in practice, fluctuations of a generic time series affected by 
intermittency, alternate intervals of very high activity to in- 
tervals of quiescence. 

Because of this lack of Gaussianity, the study of the fluc- 
tuations based on conventional spectral analysis is strongly 
limited, and the second order moment of the distribution is 
not longer the limiting order. An alternative way for char- 
acterizing the fluctuations is to investigate directly the dif- 
ferences of a fluctuating fleld over all the possible spatial 
scales and look at moments of orders higher than 2, adopt- 
ing the so-called multifractal approach (Parisi and Frisch, 
1985). A convenient statistical tool to perform this study 
is the so-called p — th order structure function (SF) de- 
flned as =< \V{x + r) - V{x)\p > and is expected to 
scale as r"'' . SFs are then computed for various orders as 
a function of all the possible scales and each order provides 
a value of the scaling exponent Sp. If observations show a 
non-linear departure from the simple Sp = p/3 (or Sp = p/A 
for the MHD case (Carbone, 1993)) this is an indication 
that intermittency is present. This method was introduced 
for the first time in space plasma studies by Burlaga (1991) 
who studied the exponents Sp of structure functions based 
on Voyager's observations of solar wind speed at 8.5 AU. 
This author found that, similarly to what is found in ordi- 
nary laboratory turbulent fluids, the exponent Sp was not 
equal to p/3, as expected in the K41 theory. This expo- 
nent was found to scale non-linearly with the order p and 
to be consistent with a variety of newer theories of inter- 
mittent turbulence, including Kolmogorov-Obukhov (1962). 
The first results obtained by Burlaga (1991) and Carbone 
et al. (1995) not only revealed the intermittent character of 
interplanetary magnetic field and velocity fluctuations but 
also showed an unexpected similarity to those obtained for 
laboratory turbulence (Anselmet et al.,1984). These results 
showed consistency between observations on scales of 1 AU 
and laboratory observations on scales of meters, suggesting 
a sort of universality of this phenomenon, which w£is inde- 
pendent on scale. 

While previous results referred to observations in the 
outer hcliosphere, Marsch and Liu (1993) firstly investi- 
gated solar wind scaling properties in the inner heliosphere. 
For the first time they provided some insights on the dif- 
ferent intermittent character of slow and fast wind, on the 
radial evolution of intermittency and on the different scal- 
ing characterizing the three components of velocity. They 
also concluded that the Alfvenic turbulence observed in fast 
streams starts from the Sun as self-similar but then, dur- 
ing the expansion, decorrelates becoming more multifractal. 
This evolution was not seen in the slow wind supporting the 
idea that turbulence in fast wind is mainly made of Alfven 
waves and convected structures (Tu and Marsch, 1993) as al- 
ready inferred by looking at the radial evolution of the level 
of cross-helicity in the solar wind (Bruno and Bavassano, 
1991). As we will see in the following, although the tools 
used in our analysis differ from those used by Marsch and 
Liu (1993) our results fully confirm their results but also add 
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some more inferences on the radial evolution of solar wind 
intermittency. 

Successively, several other papers tried to understand the 
phenomenon of intermittency in the solar wind looking for 
the best model which could fit the observations or could 
establish whether the observed scaling was closer to that 
shown by an ordinary fluid or rather by a magnetofluid as 
predicted by Kolmogorov (1941) and Kraichnan (1965), re- 
spectively. Ruzmainkin et al., (1995) studying fast wind 
data observed by Ulysses developed a model of Alfvcnic tur- 
bulence in which they reduced the spectral index of magnetic 
field fluctuations by an amount depending on the intermit- 
tency exponent. They found a close agreement with the 
expected Kraichnan scaling for a magnetofluid (3/2) and 
concluded that their results were consistent with a turbu- 
lence based on random-phased Alfven waves (Kraichnan, 
1965). 

Tu et al., (1996) re-elaborated the Tu (1988) model of 
developing turbulence including intermittency derived from 
the p-model of Meneveau and Sreenivasan (1987). They ob- 
tained a new expression for the scaling exponent that took 
into account that, for turbulence not fully developed, the 
spectral index is not defined yet. 

Carbone et al., (1995), for the first time adopted the Ex- 
tended Self-Similarity (ESS) concept (Benzi et al., 1993) 
to interplanetary data collected by Voyager and Helios, and 
looked for differences in the scaling properties between in- 
terplanetary magnetofluid and ordinary fluid turbulence ob- 
tained in laboratory. ESS is a powerful method to easily 
recover the scaling exponent of the fluctuations exploiting 
the interdependency of the structure functions of various 
orders. These authors concluded that, differences exist be- 
tween scaling exponents in ordinary (unmagnetized) fluid 
flows and hydromagnetic flows. 

Horbury and Balogh (1997) performed a comprehensive 
structure function analysis of Ulysses data and concluded 
that interplanetary magnetic fleld fluctuations are more 
Kolmogorov-like rather than Kraichnan-like. 

Vcltri and Mangeney (1999), adopting a method based on 
the discrete wavelet decomposition of the signal identified for 
the first time intermittent events. Successively, using condi- 
tioned structure-functions, they excluded any contribution 
from intermittent samples and were able to recover the scal- 
ing properties of the MHD fluctuations. In particular, the 
radial component of the velocity displayed the characteris- 
tic Kolmogorov slope Sp = p/3 while the other components 
displayed the Kraichnan slope Sp = p/4. 

All previous works dealt with the scaling exponents Sp of 
the structure functions , aiming to show that they follow 
an anomalous scaling with respect to that expected from 
K41 theory for turbulent fluids. This anomalous scaling is 
strictly related to the way Probability Distribution Func- 
tions (PDFs) of the increments change with scale. It is 
interesting to notice that if we consider fluctuations that 
follow a given scaling, say Svr =< \v{x + r) — v{x)\ >~ 
and introduce a change of scale, say r —f £r {£ > 0), we 
end up with the following transformation 5ver ^ i^dvr- The 
importance of this relation is that the statistical proper- 
ties of the left and right-hand-side members are the same 
(Frisch, 1995), i.e. PDF{5vir) = PDF{l''5vr). This means 
that if h is unique, the PDFs of the standardized variables 
5ur{x) = {v{x + r) — v{x))/ < {v{x + r) — v{x))^ >^^^ re- 
duces to a unique PDF highlighting the self-similax (frac- 
tal) nature of the fluctuations. In other words, if all the 
PDFs of standardized fluctuations Sur collapse to a unique 



PDF, fluctuations are not intermittent. Intermittency im- 
plies multifractality and, as a consequence, an entire range 
of values for h. Castaing et al. (1990) developed a model 
based on the idea of a log-normal energy cascade and showed 
that the non-Gaussian behavior of the Probability Distribu- 
tion Functions (PDF's) at small scales can be represented by 
a convolution of Gaussians whose variances are distributed 
according to a log-normal distribution whose width is repre- 
sented, for each scale r, by the parameter A^(r). This model 
has been adopted, for the first time in the solar wind context, 
by Sorriso et al., (1999) to fit the departure from a Gaussian 
distribution of the PDFs of solar wind speed and magnetic 
fleld fluctuations at small scales. As a matter of fact, Marsch 
and Tu (1994) had already shown that the PDFs closely re- 
semble a Gaussian distribution at large scales but, at smaller 
scales, their tails become more and more stretched as re- 
sult of the fact that large events have a probability to hap- 
pen larger than for a Gaussian distribution. Their results 
showed that values of A^(r) relative to magnetic field were 
higher than those relative to velocity throughout the incrtial 
range, confirming that PDF's of magnetic field fiuctuations 
are less Gaussian than those relative to wind speed fluctua- 
tions (Marsch and Tu, 1994). The same authors determined 
also the codimension of the most intermittent magnetic and 
velocity structures, suggesting that within slow wind inter- 
mittency is mainly due to compressive phenomena. More- 
over, the use of techniques recently adopted in the context 
of solar wind turbulence (Veltri and Mangeney, 1999, Bruno 
et al. 2001) based on wavelet decomposition allowed to iden- 
tify those events causing intermittency. Those events were 
identifled as either compressive phenomena like shocks or 
planar sheets like tangential discontinuities separating con- 
tiguous regions characterized by different total pressure and 
bulk velocity, possibly associated to adjacent flux-tubes. 

Lately, Padhye et al., (2001) used the Castaing approach 
to describe directly the PDFs of the fluctuations of the over- 
all interplanetary magnetic fleld components. These authors 
concluded that all the components followed a rather Gaus- 
sian statistics but they were not able to relate their results 
to those obtained by Marsch and Tu (1994) and Sorriso et 
al. (1999) who compared PDFs for different time scales. As 
a matter of fact, Padhye and co-workers referred to fluctu- 
ations respective to the mean fleld and not increments as 
it was done in the previous mentioned studies and in the 
present study. 

In this paper, we base our analysis on the concept of 
intermittency as given by Frisch (1995), following which a 
random function v{t) is said to be intermittent if the flatness 



(1) 



grows without bound as we filter out the lowest frequency 
components of our signal and consider only smaller and 
smaller scales. Thus, we will define a given time series to be 
intermittent if J" continually grows at smaller scales and, we 
will define the same time series to be more intermittent if J- 
grows faster. Moreover, if remains constant within a cer- 
tain range of scales, it will indicate that those scales are not 
intermittent but simply self-similar and, a value of JF 7^ 3 (3 
is the value expected for a Gaussian) would simply indicate 
that those scales do not have a Gaussian statistics. This is a 
simpler way than that used by Sorriso et al. (1999) to look 
at the behaviour of the flatness to infer the intermittency 
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character of the fluctuations but, what we gain in simplicity 
we loose in effectiveness to quantify the degree of intermit- 
tency and, we will only be able to evaluate whether a given 
sample is more or less intermittent than another one. In 
the following sections wc will analyze and discuss the radial 
evolution of intermittency in the inner heliospliere and on 
the ecliptic plane evaluating the behavior of .7^ as previously 
illustrated. 

2. 2. Data Analysis 

The present analysis was performed using plasma and 
magnetic field data recorded by Helios 2 during its first so- 
lar mission in 1976 when the s/c repeatedly observed the 
same corotating stream at three different heliocentric dis- 
tances on the ecliptic plane, during three consecutive so- 
lar rotations. In order to compare intermittency between 
high and low speed plasma, low speed regions ahead of each 
corotating high speed stream, were also studied. The three 
streams, named " 1" , "2" and "3" , respectively, cam be identi- 
fied in Figure 1 where the wind speed profile and the space- 
craft heliocentric distance are shown for the whole Helios 
2 primary mission to the Sun. The exact location of the 
selected intervals, lasting 2 days each, is shown by the rect- 
angles drawn on the data profile. Beginning and end of 
each time interval are shown in Table 1 where, we also show 
the average heliocentric distance, the average wind speed, 
the angle Ab-v between magnetic and velocity vectors, the 
angle Pb-r between magnetic field vector and radial direc- 
tion and the angle ^v-r between velocity vector and ra- 
dial direction. While the velocity vector is always closely 
aligned with the radial direction, magnetic field vector gen- 
erally follows the expected Archimedean spiral configuration 
although this agreement is larger during fast wind than dur- 
ing slow wind time intervals. The data set is made of 81 sec 
averages of magnetic and plasma observations recorded in 
Solar-Ecliptic reference system SE where, the X axis is ori- 
ented towards the sun, the y axis lies on the ecliptic and 
it is oriented opposite to the s/c direction of motion and, 
the 2 axis completes the right-handed reference system. 
These fast wind streams are notorious for being dominated 
by Alfvcn waves and have been widely studied since they 
offer a unique opportunity to observe the radial evolution of 
MHD turbulence within the inner heliosphere (for a rather 
complete review of existing literature related to this topic 
see Tu and Marsch, 1995). 

The aim of the present study is to investigate the behav- 
ior of magnetic field and wind velocity intermittency as a 
function of heliocentric distance and type of wind (i.e. fast 
and slow). Although intermittency refers to the statistical 
behavior of the fluctuations in the spatial domain, it can be 
estimated from measurements made in the temporal domain 
simply adopting the Taylor's frozen-in hypothesis. This as- 
sumption, which is fully acceptable within the usual condi- 
tions of strongly supersonic and super-Alfvenic solar wind, 
allows to treat, with good approximation, each fluctuation 
as an eddy and spatial r and temporal r coordinates can be 
unitually exchanged via the relation r = Vsw • t where Vsw is 
the solar wind bulk speed. In order to study intermittency 
we computed the following estimator of the flatness factor 

:f 

nr) = ^ (2) 



where r is the scale of interest and =< \V{t + r) — 
V{t)\^ > is the SF of order p of the generic function V{t). 
This definition slightly differs from that given by Frisch 
(1995) since we compute the factor for each single scale 
while Frisch calculates T using a high-pass filter whose 
cutoff frequency is repeatedly shifted towards higher and 
higher frequencies each time. However, in both cases a given 
function is considered intermittent if the factor increases 
when considering smaller and smaller scales or, equivalently, 
higher and higher frequencies. 

A vector field, like velocity and magnetic field, encom- 
passes two distinct contributions, a compressive one due to 
intensity fiuctuations that can be expressed as 

5\B{t,T)\ = \B{t + r)\-\B{t)\ (3) 

and a directional one due to changes in the vector orienta- 
tion 



SBit, r) = / ^ (i?,(f + r) - Biit)r (4) 

y i — x.y,z 

Obviously, relation 4 takes into account also compressive 
contributions and the expression 5B{t,T) > \6\B{t,T)\\ is 
always true. 

In the following we will study the flatness factor T ob- 
tained from SFs computed for both conipressive 
\S\B{t,T)\\P >t and directional C? =< {SB{t,r)y >t fluc- 
tuations. As regards this last quantity, we like to stress that 
we verified that magnetic sector changes do not apprecia- 
bly influence its value and that in this study only interval 
(72:00-74:00) contains a magnetic sector change. Compar- 
ing the radial dependence of these two quantities for fast and 
slow wind and for magnetic field and velocity will turn out to 
be useful to better interpret the radial evolution of intermit- 
tency as observed in the solar wind MHD turbulence. Our 
analysis will be based on the following definitions: 1) a given 
time series is defined intermittent if the factor !F monoton- 
ically increases moving from larger to smaller scales, 2) the 
same time series is defined more intermittent than another 
one if J- begins to increase at larger scales since, following 
Castaing et al. (1990), this implies a larger inertial range 
and, consequently, a larger number of steps along the cas- 
cade with intermittency increasing at each step, 3)if starts 
to increase at the same scale for two different time series, we 
will consider more intermittent the one for which grows 
more rapidly. Moreover, we like to remind that a Gaussian 
statistics would show values of close to 3 for all scales, in- 
dicating the self-similar character of our fiuctuations. How- 
ever, if !F fluctuates around a value somewhat different from 
3 our fluctuations are still self-similar although not Gaus- 
sian. Anyhow, in both cases these fluctuations are not con- 
sidered intermittent. 

Thus, our definition of intermittency will be limited to a 
qualitative definition rather than quantitative since the aim 
of the present work is only to compare the radial evolution of 
intermittency for different solar wind parameters and within 
different solar wind conditions. 

3. 3. Magnetic field and velocity intermittency 
vs heliocentric distance 

Values of for both scalar and vector differences for mag- 
netic field as a function of temporal scale t expressed in 
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seconds are shown in Figure 2. The factor T has been com- 
puted for slow (left colunin) and fast (right column) wind 
and for three distinct radial distances as indicated by the 
different symbols used in the plots. In addition, errors rel- 
ative to each value of T arc also shown. It is readily seen 
that magnetic field fluctuations in both slow and fast wind 
are intermittent since T increases at small scales. Values of 
T for compressive fluctuations within slow wind (panel A) 
start to increase well beyond 10* sec reaching values larger 
than 20 at the smallest scale. Moreover, there is no evi- 
dence for any radial dependence since all the curves overlap 
to each other within the error bars. On the contrary, com- 
pressive fluctuations for fast wind (panel B) show a clear 
radial dependence. As a matter of fact, the three curves 
intersect each other at large scales down to ~ 4 • 10^ sec 
but clearly separate at smaller scales indicating that inter- 
mittency increases with the radial distance from the sun. 
Moreover, since within slow wind T starts to grow at larger 
scales and reaches higher values at small scales, we can say 
that magnetic compressive fluctuations in slow wind, at least 
at 0.3 and 0.7 AU, are more intermittent than those within 
fast wind. Slow wind directional fluctuations (panel C) are 
also rather intermittent since T starts to increase around 
2-10'' sec, at frequencies slightly higher than for slow wind 
compressive fluctuations. However, these fluctuations are 
less intermittent than compressive fluctuations in the same 
type of wind since T increases more slowly at small scales. 
Moreover, there is no radial dependence. Panel D shows the 
behavior of T for directional fluctuations in fast wind. Also 
in this case as for panel B, there is a clear radial dependence 
of intermittency on the radial distance. The flatness factor 
J- remains approximately constant and rather similar for 
the three distances at large scales down to ~ 2 • 10^ sec and 
then increases more rapidly for larger heliocentric distances. 
Thus, our sample at 0.9 AU is more intermittent than that 
at 0.7 AU, which, in turn, is more intermittent than that 
at 0.3 AU. Considering that the scales at which T starts to 
increase is only around 10'^ sec and that the values reached 
at small scales are lower, these fluctuations are less intermit- 
tent than the corresponding ones within slow wind and the 
compressive ones within the same fast wind. Moreover, the 
fact that in panel D T starts to increase at much smaller 
scales than in slow wind (Panel C), is strongly indicative 
that the inertial range in this case is much less extended as 
we already know from the existing literature (see review by 
Tu and Marsch, 1995). 

Results relative to velocity fluctuations are shown in Fig- 
ure 3 in the same format adopted in the previous Figure. 
Values of T for slow speed compressive fluctuations start 
to increase around 10* sec (panel A). However, the three 
curves intersect each other various times along the whole 
range of scales showing that there is no clear radial depen- 
dence although, the smallest scale would suggest some radial 
evolution which, in addition, would be opposite to what is 
observed in fast wind. However, the large associated error 
bars do not allow us to draw any realistic conclusion. Tak- 
ing into account that the the scale at which T starts to 
increase is of the same order of that relative to magnetic 
field compressive fluctuations in slow wind but T reaches 
much lower values at small scales, we conclude that velocity 
compressive fiuctuations are less intermittent than magnetic 
compressive fluctuations in slow wind. Panel B, relative to 
velocity compressive fluctuations in fast wind, shows a clear 
radial dependence of T. The three curves start to increase 
around 10^ sec and separate at smaller scales. Then, in- 
termittency increases from 0.3 to 0.9 AU since T increases 



more rapidly for larger heliocentric distances. This result is 
similar to what we observed for magnetic compressive fluc- 
tuations in fast wind although the overall intermittency in 
this case is much reduced taking into account that T starts 
to increase at much smaller scales and roaches smaller val- 
ues. Panel C shows results relative to velocity directional 
fluctuations in slow wind. These curves, although less sta- 
ble than the corresponding curves relative to magnetic field 
(Figure 2C) , show a very similar behavior and no hint for a 
possible radial dependence. Also panel D, where we report 
values of T for velocity directional fluctuations in fast wind, 
shows results very similar to those shown for magnetic field 
in Figure 2D to the extent that these two sets of curves over- 
lap to each other, within the error. This last result, as it 
will be discussed later on in this paper, clearly derives from 
the strong contribution due to Alfvenic fluctuations popu- 
lating the fast corotating streams that we selected (Bruno 
ct al., 1985). Differently from magnetic field fiuctuations, 
velocity directional fiuctuations seem to be only slightly less 
intermittent than compressive fluctuations. 

4. 4. Intermittency in the mean field 
reference system 

Although, other authors (Marsch and Liu, 1993, Marsch 
and Tu, 1994) already addressed the study of the radial 
evolution of intermittency for magnetic field and velocity 
components, we like to provide a complete picture of this 
radial dependence adding a study performed in the mean 
field coordinate system (MF, hereafter) which, for magnetic 
field, is more appropriate than the usual RTN or SE co- 
ordinate systems. As a matter of fact, the large scale in- 
terplanetary magnetic field configuration breaks the spatial 
symmetry and introduces a preferential direction along the 
mean field. As a consequence, a natural reference system is 
the one for which one of the components, that we call B//, 
is along the mean field B outwardly oriented, and the other 
two are perpendicular to this direction. In our case we chose 
one of the two perpendicular components _B_l2 to be perpen- 
dicular to the plane identified by B and the mean solar wind 
velocity V , so that _Bj_2 = B y,V and, the remaining direc- 
tion Bxi descends from the vector product Bx2 x B, where 
the symbol ' indicates a unitary vector. In the top panel of 
Figure 4, we show values of T for the three magnetic com- 
ponents in SE reference system (i.e. Bx, By and Bz) at 
the three different heliocentric distances previously chosen 
and, in the bottom panel, we show the components in the 
MF reference system (i.e. B//, B±i and i3x2)for the same 
heliocentric distances. In the top panel increases for all 
the components as the radial distance increases. While at 
0.9 AU the three components show the same behavior, at 0.3 
and 0.7 AU the curve relative to Bx runs above the other 
two curves. Since the distance between these curves slightly 
increases at small scales we might conclude that the radial 
component is slightly more intermittent than the other two 
components. Unfortunately, this conclusion is not corrob- 
orated by the size of the errors associated to each point, 
which are quite large. The bottom panel shows results rel- 
ative to the MF reference system. At 0.3 and 0.7 AU there 
is not much difference with the situation discussed in the 
previous panel because the orientation of the two reference 
systems is not very different either, given that the magnetic 
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field is almost radially oriented (see Table 1). On the con- 
trary, moving to 0.9 AU and comparing these results with 
those obtained in the other reference system, we clearly ob- 
serve a decrease of T for the two perpendicular components 
and an increase for the parallel comijoncut. Since the throe 
curves lie on the same level at large scales and end up with 
remarkable differcut values at small scales, we conclude that 
the component parallel to the local magnetic field is more 
intermittent than the perpendicular components. Moreover, 
the two perpendicular components are less intermittent in 
the MF reference system than in SE. This means that in the 
MF reference system wo enhance on one hand the stochastic 
character of the fiuctuations perpendicular to the local field 
direction and, on the other hand, the coherent character of 
the fluctuations along the local field direction. 

In Figure 5 we show, for the slow wind, the same ele- 
ments discussed in the previous Figure. The much steeper 
behavior of these curves suggests that, generally, slow wind 
is more intermittent than fast wind. Moreover, especially at 
0.3 AU, there is a tendency for both Bx, in the top panel, 
and S//, in the bottom panel, to be steeper than the other 
components at small scales, suggesting a higher intermit- 
tency but, this tendency is not confirmed at larger heliocen- 
tric distances. As discussed in the following, the reason for 
this appreciable different behavior of Bx and Bj j might be 
due to the fact that so close to the sun the contribution due 
to Alfvenic fluctuations, mainly acting on the perpendicular 
components, is not negligible even within slow wind (Bruno 
et al., 1991). As a matter of fact, the stochastic nature of 
the fluctuations due to Alfvcn waves tends to make more 
Gaussian the PDFs of magnetic and velocity fluctuations 
perpendicular to the mean field direction. In conclusion, an 
overall view reveals that the behavior of T within slow wind 
is not very sensitive to this change of reference system. 

For sake of completeness we have rotated into the MF 
reference system also velocity fluctuations although this ref- 
erence system is not the most appropriate for this parameter 
given that the wind expands radially. The two panels of Fig- 
ure 6 show, for the three heliocentric distances, the behavior 
of T for fast wind velocity fluctuations in the SE reference 
system and in the MF reference system, respectively. We 
like to remark that in the SE reference system Vx resem- 
bles very closely the behavior of the fast wind speed shown 
in Figure 3 since the average wind velocity vector is always 
close to the radial direction. Moreover, the two perpendicu- 
lar components, Vy and Vz in SE and Vli and Vj_2 in MF, 
at 0.3 and at 0.7 AU closely recall the behavior of the corre- 
sponding magnetic components within fast wind. However, 
the presence of a large plateau in the central part of Vx and 
Vj I makes it more difficult to estimate the degree of inter- 
mittency of these components with respect to the perpen- 
dicular ones. At 0.9 AU, due to a weaker stationarity in the 
data, the situation looks even more complex and does not 
allow to estimate which component is the most intermittent 
one. T remarkably increases with distance at small scales 
for all the components in both reference systems. However, 
as expected, the rotation into the MF reference system does 
not have a laxge influence at 0.3 AU but it causes a general 
increase of T at 0.9 AU. The enhancement is such that the 
two perpendicular components have the same behavior, and 
differences with the parallel component become appreciably 
smaller. This is due to the fact that in this reference system 
the fluctuations of the components are not longer indepen- 
dent from each other as it would be in SE reference system. 

Finally, in Figure 7, we show results relative to the slow 
wind in the same format of the previous Figure. Here, the 



very confused behavior of the curves and the large associated 
errors, especially at 0.7 and 0.9 AU, suggest a rather weak 
stationarity of the data and make it difficult to compare the 
behavior of different components. A general comment that 
wo can easily make is that these curves are much steeper and 
start to increase at much larger scales than in fast wind. As a 
consequence, velocity components in slow wind are generally 
more intermittent than in fast wind. In addition, the rota- 
tion from SE to MF does not influence much our results, as 
expected. However, it is such that the behavior of the three 
velocity components at 0.3 AU looks more similar to that 
of the corresponding magnetic components (Figure 5, lower 
panel). This suggests that Alfven waves, although less rele- 
vant than in fast wind, might play a role even in this sample 
of slow wind. 

5. 5. Summary and discussion 

We studied the radial dependence of solar wind intermit- 
tency looking at magnetic fleld and velocity fluctuations be- 
tween 0.3 and 1 AU. In particular, we analyzed compressive 
and directional fluctuations for both fast and slow wind. Our 
analysis exploits the property that probability distribution 
functions of a fluctuating fleld affected by intermittency be- 
come more and more peaked at smaller and smaller scales. 
Since the peakedness of a distribution is measured by its 
flatness factor we studied the behavior of this parameter at 
different scales to estimate the degree of intermittency of 
our time series. Our general results can be summarized in 
the following points: 

1) magnetic fleld fluctuations are more intermittent than 
velocity fluctuations; 

2) compressive fluctuations are more intermittent than 
directional fluctuations; 

3) slow wind intermittency does not show radial depen- 
dence; 

4) fast wind intermittency, for both magnetic field and 

velocity, clearly increases with distance. 

5) magnetic and velocity fiuctuations have a rather Gaus- 
sian behavior at large scales, as expected, regardless of type 
of wind or heliocentric distance. 

Point 4 is particularly interesting because we found that 
both compressive and directional fluctuations become more 
intermittent with distance. As a matter of fact, if we think 
of relations 3 and 4 we eeisily realize that while intermittency 
of directional fiuctuations can be fully uncomprcssive, it is 
not possible to avoid that intermittency of compressive fiuc- 
tuations contaminates directional fluctuations. In the latter 
case, the limiting condition would be the same intermittency 
level for both kind of fluctuations. Thus, intermittency of 
directional fluctuations contains also contributions due to 
compressive fluctuations. This distinction plays an impor- 
tant role in discussing our results since the intermittency 
character of directional fluctuations reflects the contribution 
of both compressive phenomena and uncomprcssive fluctu- 
ations like Alfven waves. 

Now, there are at least two questions that we should ad- 
dress: 1) why directional fluctuations are always less inter- 
mittent than compressive fluctuations? and, 2) why only 
fast wind shows radial evolution? We can explain our ob- 
servations simply assuming that the two major ingredients 
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of interplanetary MHD fluctuations are compressive fluc- 
tuations due to a sort of underlying, coherent structure 
convected by the wind and stochastic Alfvenic fluctuations 
propagating in the wind. The coherent nature of the first 
ingredient would contribute to increase intcrmittency while 
the stochastic character of the second one would contribute 
to decrease it. If this is the case, coherent structures con- 
vected by the wind would contribute to the intermittency 
of compressive fluctuations and, at the same time, would 
also produce intermittency in directional fluctuations. How- 
ever, since directional fluctuations are greatly influenced by 
Alfvenic stochastic fluctuations, their intermittency will be 
more or less reduced depending on the amplitude of the 
Alfvcn waves with respect to the amplitude of compressive 
fluctuations. Thus, compressive fluctuations would always 
be more intermittent than directional fluctuations. 

Before addressing the second question, we like to re- 
call that several papers have already shown (sec review by 
Tu and Marsch, 1995) that slow wind Alfvenicity does not 
evolve with increasing the radial distance from the sun. As 
a matter of fact, power spectra exhibit a spectral index 
close to that of Kolmogorov and a rather good equiparti- 
tion between inward and outward modes (Tu et al., 1989). 
Thus, once the inertial range is established, the Alfvenicity 
of the fluctuations freezes to a state that, successively, is 
convected by the wind into the interplanetary space with- 
out major changes (Bavassano et al., 2001). On the other 
hand, within fast wind, turbulence is dominated by outward 
propagating Alfven waves, which strongly evolve in the inner 
holiosphoro becoming weaker and weaker during the wind 
expansion, to the extent that at 1 AU, on the ecliptic plane, 
their amplitude is much reduced and of the order of that of 
inward propagating Alfven waves. At that point, the result- 
ing Alfvenicity resembles the one already observed in the 
slow wind close to the sun (Tu and Marsch, 1990). Keeping 
this in mind, taking into account that convected structures 
experience a much slower radial evolution because they do 
not interact with each other non-Iinearly as Alfven waves 
do, considering that Alfven waves are mainly found in fast 
rather than in slow wind, it comes natural to expect that 
intermittency would radially evolve within fast rather than 
slow wind. Obviously, this would explain why directional 
fluctuations become more intermittent only within fast wind 
but would not explain why also compressive fluctuations be- 
come more intermittent within fast wind. In reality, if we 
consider that compressive events cause intermittency (Veltri 
and Mangeney, 1991, Bruno et al., 2001), we might ascribe 
this different behavior to the fact that fast wind becomes 
more and more compressive with radial distance while the 
compressive level of slow wind remains approximately the 
same, as shown by Marsch and Tu, (1990). 

Our analysis performed on the components can also help 
to understand, although partially, the topology of these con- 
vected structures. In SE reference system, fluctuations along 
the radial component are more intermittent than those per- 
pendicular to it as already found by Marsch and Liu (1993), 
although this feature, especially for the magnetic field, tends 
to vanish around 1 AU. The reason is that perpendicular 
components are more influenced by Alfvenic fluctuations 
and as a consequence their fluctuations are more stochas- 
tic and less intermittent. This effect largely reduces during 
the radial excursion mainly because the SE reference system 
is not the most appropriate one for studying magnetic field 
fluctuations. As a matter of fact, the presence of the large 
scale spiral magnetic field breaJis the spatial symmetry in- 
troducing a preferential direction parallel to the mean field. 



Consequently, we showed, that if we rotate our magnetic 
data into the mean field reference system, especially at 0.9 
AU, the intermittency of the perpendicular components de- 
creases and that of the parallel component increases. More- 
over, the two perpendicular components show a remarkable 
similar behavior as expected if they experience Alfvenic fluc- 
tuations. On the other hand, results obtained on velocity 
fluctuations suggest that a reference system with an axis 
parallel to the radial direction looks more appropriate to 
perform a similar study showing that the radial component 
seems to be the most intermittent component. 

One further observation is that generally most of the 
curves relative to velocity fluctuations came out to be less 
stable than those relative to magnetic field fluctuations and 
affected by larger errors due to a weaker stationarity with 
respect to magnetic field data. 

Finally, our results cannot establish whether magnetic 
and velocity structures causing intermittency are convected 
directly from the source regions of the solar wind or they 
are locally generated by stream-stream dynamic interaction 
or, as an alternative view would suggest (Primavera et al., 
2002), they are locally created by parametric decay insta- 
bility of large amplitude Alfven waves. Probably all these 
origins coexist at the same time and confirm that in any case 
the the radial dependence of the intermittency of interplan- 
etary fluctuations is strongly related to the turbulent radial 
evolution of their spectrum. 
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time interval 


radial distance [AU] 


< V > [km/s] 


aB-v[°] 




7v-fl[°l 


46:00-48:00 


0.90 


433 


29.4 


29.6 


2.9 


49:12-51:12 


0.88 


643 


31.4 


29.6 


2.4 


72:00-74:00 


0.69 


412 


15.2 


16.3 


2.9 


75:12-77:12 


0.65 


630 


19.4 


18.1 


1.5 


99:12-101:12 


0.34 


405 


22.9 


20.6 


2.3 


105:12-107:12 


0.29 


729 


8.2 


7.5 


1.8 



Table 1. From left to right: time interval in dd:hh, heliocentric 
distance in AU, average wind velocity in km/s, angle between 
field and velocity vectors, angle between field vector and radial 
distance and, finally, angle between velocity vector and radial 
distance. 
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Helios 2: day of 1976 



Figure 1. Solar wind speed profile during Helios 2 primary 
mission. Rectangles overlaying the plot indicate the time 
intervals selected for the analysis (see also Table 1 for more 
details on the intervals). 
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Figure 2. Flatness T versus time scale r relative to mag- 
nctie field fluctuations. The left column (panels A and C) 
refers to slow wind and the right column (panels B and D) 
refers to fast wind. The upper panels refer to compressive 
fluctuations, the lower panels refer to directional fluctuar 
tions. Vertical bars represent errors associated to each value 
of T. The three different symbols in each panel refer to dif- 
ferent heliocentric distances as reported in the legend. 
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Figure 3. Flatness J- versus time scale r relative to wind 
velocity fluctuations. In the same format of Figure 2 panels 
A and C refer to slow wind and panels B and D refer to 
fast wind. The upper panels refer to compressive fluctua- 
tions and the lower panels refer to directional fluctuations. 
Vertical bars represent errors associated to each value of J^. 
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Figure 4. Flatness T versus time scale r relative to fluc- 
tuations of the components of the interplanetary magnetic 
field observed in fast wind. The scale of the horizontal axis 
is divided in three parts all covering the same range of scales. 
Upper panel) there are three sets of curves at three different 
heliocentric distance. Within each set, different components 
are indicated by different symbols as reported in the legend. 
Components in this panel are taken in the Solar Ecliptic 
(SE) reference system where, the X axis is oriented towards 
the sun, the y axis lies on the ecliptic and it is oriented 
opposite to the s/c direction of motion and, the Z axis com- 
pletes the right handed reference system. 
Lower panel) parallel and perpendicular components in the 
Mean Field (MF) reference system. -B//, is along the main 

field B outwardly oriented, Bj_2 is perpendicular to the 
plane identified by 5 and the mean solar wind velocity V , 
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Figure 5. Flatness T versus time scale r relative to fluc- 
tuations of the components of the interplanetary magnetic 
field observed in slow wind, in the same format of Figure 4. 
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Figure 6. Flatness T versus time scale r relative to fluc- 
tuations of velocity components observed within fast wind. 
Results in the upper panel refer to components observed in 
the SE reference system while in the lower panel refer to 
solar wind components rotated into the mean field (MF) 
reference system. Format and symbols are the same used 
for Figures 4 and 5. 
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Figure 7. Flatness T versus time scale t relative to fluctu- 
ations of velocity components observed in slow wind, in the 
same format of Figure 6. 



